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Abstract 

A diffusion process is usually assumed for the phase of the order parameter of a Bose system of 
finite size. The theoretical basis is limited to the so called Bogoliubov approximation. We show 
that a suitable generalization of the Hartree-Fock-Bogoliubov approach recovers phase diffusion. 
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I. INTRODUCTION 



Bose-Einstein condensation (BEC) of weakly interacting atom systems attracted a large 
attention since the studies on liquid Helium, and a renewed interest stimulated by experi- 
mental observations of condensation of trapped alkali atoms [l|, Q, Si . 

n 

Relevant theoretical approaches have been reviewed by Griffin |4j). The condensate is 
usually considered in a coherent state whose amplitude satisfies a nonlinear Schrodinger 
equation [3, E , Q] - Quantum fluctuations determine the instability of this state, known as 
phase diffusion (PD), as analyzed in the work of Lewenstein and You [8j. This instability is 
expected in systems in restricted geometries, where a real symmetry-breaking phase cannot 
occur. In Ref. jsy, a gapless approximation, the Bogoliubov or the Hartree-Bogoliubov 
(HB), has been shown to exhibit PD, and it is also noted that a complete self-consistent 
approximation, the Hartree-Fock-Bogoliubov (HFB), has a gapped which spectrum prevents 
PD. However, it must be rejected because does not satisfy the number conservation law (J). 

In a different contest [9( , we recently described the possibility of observing BEC and quasi- 
superfluid behavior in two-mode photon systems. In that case the single mode approximation 
can be made. Thus, we obtain the zero-dimensional version of the problem of two interacting 



Bose condensates 



121 This model can be also viewed as a particular case of systems 
141 ] . We showed that in the study of the evolution of an initial 



of itinerant polaritons [131 . 
coherent state, there is, in the condensate phase, an initial time range where the effect of a 
symmetry-breaking field can be neglected, and PD should be observed. 

On the other hand, the Bogoliubov transformation amounts to a rotation in the particle 
degrees of freedom which leads to quasi-particles. It is easily shown that the rotation angle 
becomes infinite in the limit of vanishing symmetry-breaking field in the HB approximation, 



while it is kept finite in the HFB approximation [15|, [16|. Thus, quantities related to the 
condensate fluctuations, as for instance (aa) and (aW), in HB are actually diverging in the 
limit of vanishing symmetry-breaking field. As a consequence, the instability of the con- 
densate is strongly enhanced with respect to PD. In spite of the weakness of the theoretical 
background, the existence of PD is widely accepted and used to explain several experimental 
and theoretical results 



in, Q, Q . 



The aim of this work is to show that a self-consistent approximation can be introduced 
which satisfies both the requirements of finite rotation and the absence of a gap in the exci- 
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tation spectrum. This latter feature is a direct consequence of the validity of the continuity 
equation. This result can be achieved within a truncation procedure of the equations of 
motion which considers linear and bilinear quantities in the creation and annihilation oper- 
ators space. Higher order terms are taken into account by means of the factorization in the 
original picture of creation and annihilation particles operators, or alternatively, by normal 
ordering of corresponding quantities in the quasi-particle operators. 

The structure of the paper is the following. In Sec. [Til we introduce the general model 
and discuss the problem of PD in the Bogoliubov and HFB approximations. In Sec. II III we 
first establish the equations of motion truncation, and then we show how PD is mainteined 
in the extended selfconsitent theory. Then, we conclude the paper in Sec. [TV] 



II. THE MODEL 

We first introduce the problem in the context of interacting photon systems Two 
e.m. modes interacts via a tunneling term in a nonlinear medium: 

H = ujq (a\ai + a\a<2^ — w (a\a2 + a\ai \ + g (a\ai + a\a^ (1) 

Through a canonical transformation H reduces to 

H = (co - w) a ] a + (u + w) b ] b + g (a ] a + b ] b) 2 (2) 

where a — (ai + (12) /v2 and b = (a% — 02) / \/2. If the value of the hopping constant w 
exceeds the mode energy ujq, the the vacuum for the mode a is unstable. Under these 
conditions, by applying the Bogoliubov approximation to a coherent state for the mode a, 
a phase diffusion phenomenon appears [8j|. The interaction between a and b can be seen 
as a particular case of interaction between condensate and quasi-particles in Bose particle 
systems in the limit where only one quasi-particle mode is considered. Taking into account 
that the ground state of the isolate b system is the vacuum, we concentrate ourselves on 
the the parts of H concerning only the mode a: H a = u' a^a + ga) 2 a 2 — A (a + aH , where 
lo'q = u — w + g, and where we are considering explicitly a symmetry-breaking field. This 
field does not represent a purely mathematical tool introduced to describe the emergence 
of a superfluid phase, as usual in boson particles systems, but it can be physically realized 
through the interaction of the e.m. modes with a "classical" non-fluctuating electron current 
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[201 ] . An exact solution for the evolution of a coherent state is not known in the presence 
of the symmetry-breaking field, and different approximations can be made. A constraint 
which should be verified by any approach is represented by the continuity equation 

= _, A (a - at) (3) 
dt v ' w 

We observe first that in the limit of vanishing nonlinearity the ground state reduces to a 
coherent state, that is to the vacuum of a Hamiltonian where the degrees of freedom have 
been translated. Then, having in mind a weak-coupling theory, we perform the translation 
a — > a + v. The new Hamiltonian is 

H = E + Hi + H 2 + H 3 + H 4 , where 



Eq 


= uj' u 2 + gv 4 — 2Xu, 




(4) 


Hi 


= (a f + a) (uj'qV + 2#z/ 3 - 


A), 


(5) 


H 2 


= uj'qo)cl + gu 2 (4n a + a} 2 


+ a 2 ) , 


(6) 


H 3 


= 2gu (a) 2 a + a^a 2 ) , 




(7) 


H 4 


t2 2 

= ga' a . 




(8) 



In this new representation the continuity equation reads as 

da^a ( da da' 



( da da^\ , + x 



dt 

The occurrence of BEC implies gv 2 finite even for small g. 

The Bogoliubov approximation amounts to take into account terms of order gv 2 and to 
neglect terms of order ^fg and g. In this limit we can disregard H% and H^. The condensate 
amplitude v can be fixed by minimizing E Q : 

u' Q v + 2gv i - A = 0. (9) 

A finite solution for the condensate amplitude u, in the limit of vanishing A, is obtained only 
for oj'q < 0. We note that this conditions makes H\ vanishing. Due to this choice, 

H = -a ] a + gv 2 (2a j a + a j2 + a 2 ) . (10) 

For \jv gv 2 the ground state is in an eigenstate of the quadrature x = {a) + a) /y/2, 
while the other quadrature p = i (a^ — a) / y/2 will have infinite fluctuations. As far as the 



evolution of an initial coherent state is concerned, we obtain PD. Indeed, the coherent state 
of amplitude a evolves as 



l«)t= < 



iy/2ap{t) 



|0) = exp \iV2a [p(0) - AigvHx (0)]1 |0) . 



(11) 



The corresponding wave function in the x representation is then 



$ a (x, *) = 4 ex P 



ax — — (l — Aigu 2 t) 



(12) 



On the other hand, we can characterize PD in terms of the following average quantities 



(a\ x (t) \a) = xo 
(a\p (t) \a) 



(a\ p 2 (t) \a) — (a | p (t) \a) 



-Agv x t 
- + 8g 2 u 4 x 2 t 2 



(13) 
(14) 
(15) 



The HFB approximation amounts to treat H 3 and in a mean-field approxi- 



mation. If one admits that hig. 



his 


ier 


15 




16 



ier order terms can correct H\ and H2 and lead 



to a squeezed ground state |15|, ll6] which can be defined as the vacuum of 7 = 
(cosh#a — sinh#at), H 3 becomes Igv (2 sinh 2 9 + cosh 9 sinh 9) (a + a^), while H4 is now 
g [4sinh 2 9a) a + cosh^ sinh 9 [a) 2 + a 2 )] . These assumptions being made, the angle 9 is 
given by the self-consistent equation 

2g {y 2 + sinh 9 cosh 9) 



tanh 29 



(16) 



2g(v 2 - sinh 9 cosh 9) + £ 

which admits a finite value. While the Bogoliubov approximation satisfies the continuity 
equation at least in the average, this is not true for HFB, because dx/dt is proportional to 
g sinh 9 cosh 9. 



III. THE EXTENDED HFB SOLUTION 



Now we introduce a mean field approach going beyond the simple Bogoliubov approxima- 
tion without introducing a gap in the excitation spectrum. We discuss the problem starting 
from the equations of motion. First, we assume the existence of a finite angle 9 and perform 
the Bogoliubov transformation. The Hamiltonian is 
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Hi = A (7 + 7 f ) 

H 2 = S 7 f 7 + 7; Ei (tV + 77) 



#3 = 201/ 

#4 = 4 

where 



(7 3 + 7 t3 ) %; sinh 20 + (7+77 + 7V7) e e ('cosh 20 + 



(17) 
(18) 

(19) 



^cosh 2 20 + ^ sinh 2 20^ 7 t2 7 2 + J sinh 2 20 ( 7 t4 + 7 4 ) + cosh 29 sinh 20 ( 7 t 7 3 + 7%)) ) 



A = z/e 6 



(cud + 2 ^ 2 ) + 2 £ sinh e (e 9 + sinh 0) 

£ = A + 3 [cosh 2 20 + sinh 2 20 + cosh 20 (2 sinh 2 0-1)] 
A = (c^ + 4#z/ 2 )cosh20 + 2#z/ 2 sinh20 
E 1 = Ai + 2# sinh 20 ^cosh 20 + sinh 2 9 - i 
Ai = (c^ + 4#z/ 2 ) sinh 20 + 2#z/ 2 cosh 20 



(21) 
(22) 

(23) 
(24) 



Then, we build a set equation of motions by limiting ourselves to consider only linear and 
bilinear terms and neglecting higher order operators: we choose the operators 7, 7', 7 2 , 7^, 
and 7^7 as independent variables. The set of coupled equations reads 

i^y- = A + £07 + Eirf + 3gue 6 sinh 20 7 t2 + gue 6 (2 cosh 20 + sinh 20) ( 7 2 + 27M (25) 
at 

i^f- = _A _ ^ o7 t _ £ l7 _ 3 ffi/e e sin h 2 fl 7 2 - 3^ (2 cosh 20 + sinh 20) ( 7 t2 + 27+7) (26) 



'~dt 



Ei + 2 [A + g^e e (2 cosh 20 + sinh 20)] 7 + 6£fz/e e sinh 29^ 

1 



+2 



£ + g ( cosh 2 20 + - sinh 2 20 



7 2 + 3# sinh 2 207 



12 



+2 (Ei + 3g cosh 20 sinh 20) 7+7 



(27) 



rf 7 t2 
%—— 
dt 



-Ei - 6gue e sinh 207 - 2 [A + gue e (2 cosh 20 + sinh 20)] 7 1 " 



-3#sinh 2 207 2 -2 
-2 (Ei + 3g cosh 20 sinh 20) 7+7 
.g?7 7 



E + g ( cosh 2 20 + - sinh 2 20 



7 



12 



= A( 7 t - 7 )+£i ( 7 t2 -7 2 ) 



(28) 
(29) 
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Note that the truncation is performed after normal ordering of higher order terms. Then, 
the equations for the bilinear operators are not equivalent to those of the linear ones, i.e., 
for example, 7 2 is not simply the square of 7. 

The natural choice for the parameters v and 9 cancels all the constants appearing in the 
above equations which should give rise to instability also for A 7^ 0. This choice corresponds 
to the constraints A = and E\ = 0. Due to these conditions 

E = - 2g sinh 29^j e~ 2d . (30) 

An important consequence of these constraints is the conservation of the quasi-particles 
number n 1 = 7^7: 

^ = (31) 
dt v ; 

Furthermore, the value of 9 is exactly that of Eq. ( |T6l) . 

The important feature of the previous approximation is that the continuity equation is 
valid independently on the satisfaction of A = and E\ = 0. Once the constraints are taken 
into account, the continuity equation implies 

.sinh 29 d , 9 +9 n .ad 



i- 



, f ( 7 2 + 7 t2 ) + iue 9 ^ (7 + 7 f ) = Xe- e (7 - 7 f ) ■ (32) 



The proof requires some tedious algebra which will be reported in appendix. It is worth to 
note that for A = a new constant of motion arises which is given by (sinh 29/2) (y 2 + 7^) + 
ve e (7 + 7 f ) . 

In the Laplace space, the equations of motion can be summarized as follows: 

4 

uTi (u) = - 22 fok^k ~ (t = 0) (33) 
fe=i 

where = 7, T 2 = 7^, T 3 = 7 2 , and T 4 = 7"^, while the continuity equation implies 

sinh 29 a , , , . A 9fl . . 

— e- e 03i + 04i + 0ii + 02i = -e- 2e 34 

Iv V 

sinh 29 a , , , , , A 9fl ,„^ s 

— e~ e 032 + 042 + 012 + 022 = ^ 35 

2v v 

+013 + 023 = (36) 

sinh 29 a , , , , , „ 

-e~ 6 (034 + 044) + 014 + 024 = (37) 



2v 
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From these equations we see that it is possible to substitute the first row of the determinant 



where D ik is the minor associated to the matrix element uj5ik + <pik- It is immediately verified 
the existence of a pole in the origin (u — 0) for A — > 0. This result is expected on the basis 
of general arguments (Goldstone theorem). It is worth to note that D ik (cu) are finite in the 
present approximation. 

The resiliency of PD is directly related with the presence of the pole in the origin. Indeed, 
by considering the evolution of p (t) we find 



where fi are some nonvanishing functions of g, is, and 6. The presence of the term propor- 
tional to t is responsible for PD. Due to this term, the variance of p on a coherent state \a) 
grows as t 2 , as in the Bogoliubov case. 

IV. CONCLUSIONS 

In this paper we discussed the problem of phase diffusion in finite-size systems which 
exhibit Bose-Einstein condensation. We revised the Bogoliubov and the Hartree-Fock- 
Bogoliubov approximations. While the first one satisfies at least weakly the constraint 
of the continuity equation, the latter violates this constraint. Then, we extended the self- 
consistent approach introducing a new approximation based on the method of equations 
of motion which satisfies the continuity equation. This approximation is developed in the 
space of Bogoliubov quasi-particles. The merit of the present method is to reconcile the 
existence of a gapless spectrum with the absence of divergencies in the vanishing symmetry- 
breaking limit. At the present, while the extension to systems with many degrees of freedom 
seems to be straightforward, the validity of PD in higher order truncation schemes is under 
investigation. 
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of coefficients 4>ik with the {(A/z/)e 
space reads 



(A/V)e 26l ,0,0}. The evolution matrix in the u- 




(38) 



p (t) = {hx + f 2 x 2 + hp 2 ) t + Up + f 5 (xp + px) 



(39) 
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APPENDIX 



In this appendix we show that the continuity equation is satisfied within the approxima- 
tion performed in this paper. The two independent conditions to be fulfilled are 

sinh26 l a , . A 



2v 
sinh 29 



2v 



?32 + 042) + 012 + 022 
^33 + 043) + 013 + 023 



-— e 
v 



-26 



0. 



By writing explicitly the coefficients 
sinh 29 



A 



2v 



e~ [2gve e (2 cosh 29 + sinh 29) - Qgise 9 sinh 29] + E = -e 



-2(9 



(A-l) 
(A-2) 



(A-3) 



sinh 29 



2v 



-e-° I 2 



E + g ( cosh 2 29 + - sinh 2 29 



-3g sinh 2 29 

+gve ti (2 cosh 29 + sinh 29) - 3gve e sinh 29 = 0. (A-4) 



The first one reduces simply to 

2gsmh29e~ w + E = -< 

v 

which gives exactly the value of E written in Eq. (!30l) . The second one reads 

sinh 29 (E + g) + 2gv 2 = 0. 

To show that this equality holds we use the condition A = putting it in Ei 
lead to 



(A-5) 



2g^ + g sinh 29e~ Ae + - sinh 29e 

v 



A 



-20 



0. 



Due to this condition, Eq. 
A 



becomes 



sinh 29 I -e' 20 - 2g sinh 29e~' M + g I - - sinh 20e~' M - g sinh 29e 



-26 



V 



A 



-26* 



-10 



V 



(A-6) 



0, which 



(A-7) 



(A-8) 



Now, it is simple to show that the left hand side vanishes. 
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